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Abstract — We focus on the detection of communities in 
multi-scale networks, namely networks made of different 
levels of organization and in which modules exist at 
different scales. It is first shown that methods based on 
modularity are not appropriate to uncover modules in 
empirical networks, mainly because modularity optimiza- 
tion has an intrinsic bias towards partitions having a 
characteristic number of modules which might not be 
compatible with the modular organization of the system. 
We argue for the use of more flexible quality functions 
incorporating a resolution parameter that allows us to 
reveal the natural scales of the system. Different types 
of multi-resolution quality functions are described and 
unified by looking at the partitioning problem from a 
dynamical viewpoint. Finally, significant values of the 
resolution parameter are selected by using complementary 
measures of robustness of the uncovered partitions. The 
methods are illustrated on a benchmark and an empirical 
network. 

Index Terms — community detection, complex networks, 
modularity, multi-scale. 

I. Introduction 

Many systems of current scientific interest are made 
of elements in interaction and can be represented as 
networks. Important examples include the Internet, tele- 
phone networks, collaboration networks, airline routes, 
but also a wide range of biological networks, such as 
food-webs, metabolic networks and protein interaction 
networks. The mathematical and empirical study of 
networks has emerged in the last decade as one of 
the fundamental building blocks in the wider study of 
complex systems CO, HI , Q- One of the main reasons 
for this success is the possibility to analyze systems 
of a very different nature within a single framework. 
This approach allows to uncover similarities between 
the structures of various complex systems, which can 
reveal the existence of generic organization principles. 
A good example is the omnipresent multi-scale modular 
organization of complex networks, namely the fact that 
they are made of modules at different scales (see Fig. 1). 




Fig. 1. Sketch of a multi-scale weighted network. The width of the 
links is proportional to their weight. This network is clearly made 
of modules at different scales: 8 single nodes, 4 pairs of strongly 
connected nodes, 2 groups of 4 nodes and the system as whole. This 
multi-scale network is hierarchical as modules at one level are nested 
into modules at the next level, but this is not necessarily the case, 
i.e., some multi-scale networks are not hierarchical. 



Modules, also called communities, are defined as sub- 
networks that are locally dense even though the network 
as a whole is sparse |4]. The presence of modules at 
different scales is known to confer a crucial evolutionary 
advantage and to accelerate the emergence of complex 
systems by providing stable intermediate building blocks 

El. 

The capacity to collect large data-sets of relational 
data has radically changed the way networks are con- 
sidered and has led to the development of statistical 
methods for the description of their multi-scale topology 
and the detection of significant connectivity patterns. 
A powerful set of methods consists in uncovering the 
modules present in the network |6), Q. This identifi- 
cation has the advantage of providing a coarse-grained 
representation of the system, thereby allowing to sketch 
its organization and to identify sets of nodes that are 
likely to have hidden functions or properties in common. 
Most community detection methods find a partition of 
the nodes into communities, where most of the links 
are concentrated within the communities. Each node is 
assigned to one and only one community, i.e., partitions 
are not compatible with overlapping communities (H, 
@ . At the heart of most partitioning methods, there is a 



mathematical definition for what is thought to be a good 
partition. Once this quality function has been defined, 
different types of heuristics can be used in order to find, 
approximatively, its optimal partition, i.e., to find the 
partition having the highest value of the quality function. 

In this article, we first describe different multi- 
resolution quality functions, namely quantities incor- 
porating a resolution parameter allowing to tune the 
characteristic size of modules in the optimal partition. 
We show that these quantities are linearised versions of 
a quality function called stability iPTOl , which is based on 
the exploration of the network by a random walker at dif- 
ferent time scales. Finally, we focus on the optimization 
of these quality functions and on the important issue of 
detecting significant values of the resolution parameter 
in practical applications. The methods are successfully 
tested on a benchmark and on a real-world network. 

II. Modularity and its limitations 

Let A be the adjacency matrix of a weighted, undi- 
rected network. A is therefore symmetric and Aij is the 
weight of the link between i and j. The strength of node 
i is defined as k{ = J2j Aij\ m = Mjl^ is the total 
weight in the network. If the network is unweighted, ki 
and m are the degree of node i and the total number 
of links respectively. The quality of the partition of a 
network is a function of the adjacency matrix A and 
of the partition V of the nodes into communities. The 
widely-used modularity ifTTTl of a partition V measures if 
links are more abundant within communities than would 
be expected on the basis of chance 

Q = (fraction of links within communities) 

— (expected fraction of such links) (1) 



and reads 
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(2) 



where i,j G C is a summation over pairs of nodes 
i and j belonging to the same community C of V 
and therefore counts intra-community links. The null 
hypothesis is an extra ingredient in the definition and 
is incorporated in the matrix P^. P^ is the expected 
weight of a link between nodes i and j over an ensem- 
ble of random networks with certain constraints. These 
constraints correspond to known information about the 
network organization, i.e., its total number of links and 
nodes, which has to be taken into account when assessing 



the relevance of an observed topological feature. Two 
standard choices for the corresponding null models are 

Pij = (k) 2 /2m, then Q = Q unif (3) 

where (k) = 2m /N is the average strength and the only 
constraint is thus the total weight in the network, and 



Pij = kikj/2m, then Q = Q 



conf ■ 



(4) 



where randomized networks now preserve the strength 
of each node. The latter null model is usually preferred 
because it takes into account the degree heterogeneity 
of the network lfl2l . More complicated null models can 
in principle be constructed in order to preserve other 
properties of the network under consideration lfT3l . fl4l . 
It is interesting to note that Q un if and Q C onf are 



ij 
ij 



naturally related to the combinatorial Laplacian L 
Ay — kidij and the (normalized) Laplacian L. 
Aij/kj — 5ij respectiveljQ, and, more generally, to the 
dynamics induced by these operators (see section Ull-Cb - 
For Qconf > this relation is particularly clear after express- 
ing modularity in terms of the (right) eigenvectors v a of 
Lij, i.e., v a satisfy J2j LijV a j = \ a v a ,i- Without loss of 
generality, we assume that Ai > A2 > ... > X a > ... > 
Xn- The dominant eigenvector v± of eigenvalue Ai = 
is given by v\-i = ki/2m and is unique if the network 
is connected. By using a spectral decomposition of L^, 
one finds flOl 
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N 



conf 
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(5) 



where the contribution of the dominant eigenvector v\ 
and the null model have cancelled each other out. 

The optimization of modularity has the advantage of 
being performed without a priori specifying the number 
of modules nor their size. This procedure has been shown 
to produce useful and relevant partitions in a number 
of systems |fl"5l . Unfortunately, it has also been shown 
that modularity suffers from several limitations, partly 
because modularity optimization produces one single 
partition, which is not satisfactory when dealing with 
multi-scale systems. Related to this issue, there is the 
so-called resolution limit of modularity |fT6*1 , namely 
the fact that modularity is blind to modules smaller 
than a certain scale. This point originates from the 
bias of modularity towards modules having a certain 
scale which might not be compatible with the system 



Strictly speaking, the normalized Laplacian of a network is L i 
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but L and L are equivalent by similarity as 



architecture. This incompatibility also makes modularity 
inefficient in practical contexts as it may lead to a high 
degeneracy of its landscape lTl7l . i.e., the existence of 
several distinct partitions having a modularity close to 
the optimum, which implies that approximate solutions 
of the optimization problem are very dissimilar and that 
a partition derived from modularity optimization has to 
be considered with caution. 

III. Multi-scale methods 

A. Local maxima of modularity 

Different methods have been proposed to go beyond 
modularity optimization. A first set of methods looks 
for local maxima of the modularity landscape in order 
to uncover partitions at different scales lfl"8l . A good 
example is the so-called Louvain method, which is a 
greedy method taking advantage of the hierarchical or- 
ganization of complex networks in order to facilitate the 
optimization of modularity lfl9l . This heuristic performs 
the optimization in a multi-scale way: by comparing the 
communities first of adjacent nodes, then of adjacent 
groups of nodes found in the first round, etc. It has been 
shown in several examples that modularity estimated by 
this method is close to the optimal value obtained with 
slower methods, but also that intermediate partitions are 
meaningful and correspond to communities at intermedi- 
ate resolutions EDI . 11211 . This approach has the advan- 
tage of being fast, but it lacks theoretical foundations 
and is not able to uncover coarser partitions than those 
obtained by modularity optimization. Moreover, it may 
produce hierarchies even when the system is single-scale 
or, worse, completely random (see [20] for a discussion 
of how to deal with this issue). 

B. Multi-resolution quality functions 

Another class of methods is based on multi-scale 
quality functions. These quality functions incorporate a 
resolution parameter allowing to tune the characteristic 
size of the modules in the optimal partition and aim at 
uncovering modules at the true scale of organization of 
a network, i.e., not at a scale imposed by modularity 
optimization. The two most popular multi-scale qual- 
ity functions are ad-hoc, parametric generalizations of 
modularity. A first quantity is the parametric modularity 
introduced by Reichardt and Bornholdt (22), |23) 

^7 = ^- E E \ A n - -r p ij\ > 

which is usually defined for the configuration null model 
Pij = kikj/2m and mainly consists in changing the 



effective size of the system m c g = 771/7. The optimiza- 
tion of Q y leads to larger and larger communities in the 
optimal partition when 7 is decreased. This approach 
makes use of the size dependence of modularity: because 
of the factor l/2m in the null model, modularity depends 
on the total size of the network and not only on its local 
properties^. Decreasing 77i c fj (increasing 7) increases the 
expected weight yPy of a link between i and j, which 
makes it less advantageous to assign i and j to the same 
community (because Aij — jPij decreases). 

An alternative approach proposed by Arenas et al. ll24ll 
keeps modularity unchanged but modifies the network by 
adding self-loops to the original network. This approach 
therefore consists in optimizing 



Q, 



Q(Aij + rlij) 



(7) 



As expected, increasing r has a tendency to decrease the 
size of the communities and the optimal partition of Qoo 
is made of single nodes. Even if increasing 7 and r has, 
qualitatively, the same effect on the characteristic size 
of the communities, one should keep in mind that Q 7 
and Q r are in general optimized by different partitions, 
except if the network is regular and the resolution param- 
eters verify 7 = 1 + r/(k). It is also interesting to note 
that the quality function {/J was first proposed in order 
to preserve the eigenvectors of the adjacency matrix, as 
the eigenvectors of + rl^ and A^ are obviously 
the same. From a partitioning viewpoint, however, the 
eigenvectors of Aij do not matter as much as the eigen- 
vectors of the combinatorial Laplacian L^p ll25l and the 
normalized Laplacian Lij [26]. Moreover, modularity is 
related to the eigenvectors of the Laplacian and not of 
the adjacency matrix, see ©. These observation suggest 
to adapt the unfitting quality function {/J and to optimize 
the modularity of a modified adjacency matrix preserving 
the eigenvectors of L^. This can readily be done by 
adding strength-dependent self-loops to the nodes 



A 



13 



and by optimizing the quality function 

W 



Q r = Q(A ij + r 7 ^rS ij ) 



(8) 



(9) 



In a nutshell, this size dependence originates from a choice of 
null model where each pair of nodes i and j can be connected, given 
a certain number of available links m in the system, whatever the 
distance between i and j in the network. Local null models where 
pairs of nodes are randomly connected only within a finite radius of 
interaction are expected to suppress this effect. 



This quality function is equivalent, up to a linear trans- 
formation, to Q 1 for any network, i.e., not only for 
regular networks, with 7 = 1 + r/(k), thereby providing 
two alternative interpretations to resolutions parameters. 

C. Stability 

The multi-resolution quality functions defined in the 
previous section have been successfully tested on multi- 
scale benchmark and empirical networks |[23l , |[27ll , [28 1. 
They have the further advantage of being mathemati- 
cally very similar to modularity and of being optimized 
by modularity optimization algorithms with minimum 
code development. Unfortunately, the introduction of 
a resolution parameter, 7 or r, feels like a trick and 
lacks theoretical ground. In order to define a resolution 
parameter in a more satisfying way and, as we will see, 
to provide a more solid foundation to Q 7 and Q r , we 
look at communities from a different angle, not from 
a combinatorial point of view, where intra-community 
links are counted as in (f2l), but from a dynamical point 
of view. 

Our starting point is the following: a flow taking place 
on a network is expected to be trapped for long times 
in good communities before being able to escape [29], 
[10]. This argument suggests to measure the quality of a 
partition in terms of the persistence of flows taking place 
on the network ifTOll , |[30l . Without loss of generality, we 
describe a stationary Markov process M. as a random 
walk process. Under the condition that A4 is ergodic, 
i.e., any initial configuration asymptotically reaches the 
unique stationary solution, stability is defined as 

Rm{^) = (probability for a random walker to be in the 
same community initially and at time t) 
— (probability for two independent random 
walkers to be in the same community) (10) 

when the system is at equilibrium. 

In order to clarify this general concept, let us focus 
on a generic Markov process |[30l , namely a continuous- 
time random walk where waiting times are indepen- 
dent, identical Poisson processes. The density of random 
walkers on node i at time t, denoted by Pi(t), evolves 
according to the rate equation 



Pi 



Pj - Pi, 



(ii) 



where 



An 



~Y — uij = Lij is the Laplacian operator described 
above. In this unbiased process, a walker located at j fol- 
lows a link going to i with a probability proportional to 



Aij. If the network is connected, the stationary solution 
is unique and given by the dominant eigenvector of L, 
namely p* = ki/2m. By definition ([TOl . the stability of 
a partition associated to the Markov process (fTTb is 



(12) 



R NL{t) = yy ( e tL) — ^ Ki_Kj_ 

w 4^.. JV 'v 2m 2m 2m 

where NL stands for Normalized Laplacian. This ex- 
pression clearly shows that stability depends on time. 
The quality of a partition is thus measured differently 
at different time scales and is, in general, optimized by 
different partitions when time is tuned, thereby leading 
to a sequence of optimal partitions. 

By looking at limiting values of t, one can show that 
time acts as a resolution parameter ifTOl . fl3Ql . As time 
grows, the characteristic size of the communities is thus 
adjusted to reveal the possible multi-scale organization 
of the system. In the limit t — > 0, keeping linear terms 
in t in the expansion of -Rnl(^) leads to 

RNh(t) « (l-t)-%L(0)+tQconf = Qnl(*), (13) 

which is equivalent up to a linear transformation to Q 7 
and Q r when Pij = kikj/2m (with t = I/7, t = 
(k)/(r + (k))). These multi-resolution quality functions 
can therefore be seen as a simple linear approximation 
of R^i,(t), which provides a physical interpretation to 
the resolution parameter r and 7, i.e., the inverse of the 
time used to explore the network. It is also interesting to 
note that the configuration null model naturally emerges 
from the definition of stability and from the dynamics 
(fTTb . Interestingly, other null models, including the uni- 
form null model, are associated to other random walk 
processes 11301 . In the limit t — > 00, making use of the 
spectral decomposition of L, stability simplifies as 

flNL(t) « ^e* A2 E E '•2:«''2:,- (14) 

171 c i,jeC 
where it is assumed that the second dominant eigenvalue 
A2 of L is not degenerate and v% is its corresponding 
(right) eigenvector. -Rnl(0 is therefore maximized by a 
partition into two communities in accordance with the 
normalized Fiedler eigenvector |[26l . 

.Rnl(*) differs from modularity in several ways OOl . 
However, one can show that -Rnl(£) is always equal 
to the modularity Q con f of a time-dependent weighted 
network whose adjacency matrix is X{j(t) = (e tL ^j kf. 

RNL(t) = QcanftM*))- ( 15 ) 

This new network is symmetric if the original network is 
symmetric and the weight on its links corresponds to the 



number of walkers going from j and i in time t, when the 
system is at equilibrium. By construction, X(t) is more 
and more extended when t is increased. The optimization 
of its modularity is therefore expected to uncover larger 
communities. After noting that the Laplacian of Xy and 
the Laplacian of A+j have the same eigenvectors, (fT5T ) 
emphasizes that i?NL(t) naturally fits the arguments used 
to define Q r , see ©. 

IV. Optimization, robustness and selection of 

SIGNIFICANT SCALES 

Let us now discuss the practical side of this work, 
namely the detection of multi-scale communities in large 
empirical networks. In what follows, we will focus on 
the optimization of Qnl(^) (and equivalently of Q 7 and 
Q r ), while keeping in mind that the optimization of the 
full stability i?NL (*) can be performed by using spectral 
or greedy methods ll30l . Depending on the size of the 
network under consideration, generalizations of different 
modularity optimization heuristics can be used in order 
to optimize Qnl(*)> such as Simulated Annealing, Spec- 
tral Methods or Greedy Methods for small (N < 10 2 ), 
intermediate (N ~ 10 3 ) and large (N > 10 4 ) sparse 
networks respectively. In the following, we perform the 
optimization of Qnl(*) by using a generalization of 
the Louvain method |[T9l mentioned aboveH. One should 
stress that the outcome of the algorithm is deterministic, 
except in the initial ordering (labeling) of the nodes. This 
implies that different optimal partitions (local maxima of 
Qnl(*)) can in principle be uncovered when the initial 
ordering is changed. 

Partitions at different values of t are found inde- 
pendently by optimizing Qnl(*)> thereby producing a 
sequence of partitions that are optimal at different scales. 
However, one expects that only a small number of these 
partitions are significant, which raises another question: 
how can one select the most significant partitions, or 
equivalently the most significant scales of description 
of the network? It is ironical to note that we are thus 
confronted with a problem similar to the one that initially 
led to the definition of modularity. Modularity was 
indeed first proposed to find the best partition in a nested 
hierarchy of possible community divisions ifTTI . As we 
have argued before and will show below on an example, 
modularity does not appropriately detect important scales 
of description. 

In order to address this problem, it has recently been 
proposed to look for robust partitions, where robustness 

3 Codes are available on http://www.lambiotte.be. 



has been defined differently by different authors. This 
approach formalizes the intuitive idea that a significant 
partition should not be altered by small modifications. 
A standard measure to compare two partitions V\ and 
V2 is the so-called normalized variation of information 
V(Vi,V2) il3ll . which is a number between and 1 and 
is equal to only when the partitions are identical. Three 
types of modifications have been proposed: 

Modifying the network by reshuffling a fraction of 
the links 11321 or randomly perturbing the weight of the 
links [331 - In the following, we implement the second 
approach by randomly adding ±10% to the weight of the 
links. In practice, we optimize <3nl(£) f° r K different 
realizations of the perturbed network for each value of t, 
by always using the same node ordering. The robustness 
of partitions at time t is given by 

(vu(t) = 2 _ J2 E V(V k (t),V k >m 

^ > k=l k'=k+l 

(16) 

where Vk(t) is the optimal partition of the k th realization 
of the perturbed network at time t. In this approach, 
a scale is significant when a small modification of the 
network does not alter too much the partition found by 
the optimization algorithm. 

Modifying the optimization algorithm by taking 
advantage of the dependence of the algorithm on its 
initial condition, i.e., the node ordering |[27l . EDI . To 
do so, we optimize Qnl(*) of the original network T 
times by attributing a different, random ordering to the 
nodes. Robustness at time t is 

T T 

(vU so (t) = E v(r t (t),vAt)), 

^ > i=i i'=i+i 

(17) 

where Vi(t) is now the optimal partition when using 
the i th random ordering of the nodes at time t. In this 
approach, robustness measures the size of the basin of 
attraction.] of the optimal partitions. 

Modifying the quality function by tuning the reso- 
lution parameter [28], |[30l . To do so, we perform one 
optimization of Qnl(^) for each t, while keeping the 
node ordering fixed throughout the different values of t. 

4 The deterministic optimization process can be seen as a trajectory 
in the space of partitions toward a (ideally global) maximal partition. 
The optimization always starts from the finest partition, where each 
node belongs to its own community, but its next steps depend on 
the ordering on the nodes. The basin of attraction of an uncovered 
partition is the set of initial node orderings that lead to it. 



The robustness of partitions at time t is 



(V) QF (t) 



(18) 



where V(t) is the optimal partition at time t. In this 
approach, robustness corresponds to the persistence of 
an optimal partition over long periods of time, and to the 
fact that optimal partitions are weakly altered by tuning 
t. 

In each case, robustness is related to the ruggedness 
of the quality function landscape. Lack of robustness 
corresponds to high degeneracy, namely to the existence 
of incompatible partitions that are local maxima of 
RNL,(t) with a value close to the global maximum. 
Significant partitions are uncovered by identifying values 
of the resolution parameter where these measures of 
robustness are significantly low. 

V. Tests of the method 

In this section we test these ideas by focusing on a 
computer-generated network and a real-world network 
for which the community structure is already known. In 
each case we find that the method reliably detects the 
known community structure and reveals the important 
scales of description. 

A. Hierarchical benchmark network 

This randomly-generated network is made of 640 
nodes with 3 known hierarchical levels: small modules 
of 10 nodes nested in medium-size modules of 40 
nodes themselves nested in large modules of 160 nodes 
[18]. The expected number of links across modules and 
therefore the sharpness of the modules is tuned by a 
single parameter p, p = 1.0 in this example. We focus 
on one single realization of this random network. In 
order to evaluate {V) nct (t), (V)^ go (t) and (V}Q F (t), 
we use K = 10, T = 10 and A = 5. As shown in 
Fig. 2, the method clearly uncovers the correct scales of 
description and only those scales. The natural partitions 
into 4, 16 and 64 modules respectively are characterized 
by regions of t where (V)(t) -C 1. Moreover, these 
regimes are clearly separated by peaks of (V)(t), i.e., 
values of the resolution parameter where the algorithm 
finds conflicting partitions. It is interesting to note that 
(V)net(t), (V)aigo(*) and (V)Qp(t) have similar but non- 
identical patterns, which suggests to combine the use of 
different notions of robustness in order to improve the 
detection of significant partitions. Partitions uncovered 
by modularity (t = 1) optimization have, on average, 
15.9 ± 0.3 modules when measured over 100 optimiza- 
tions with random orderings. 




time 



Fig. 2. Analysis of a hierarchical network with N = 640 nodes 
and 3 known hierarchical levels (natural partitions into 4, 16 and 
64 modules). In the upper figure, we plot the average number of 
modules as a function of the resolution parameter t when performing 
100 optimizations with random orderings (associated to the definition 
of (^)aigo(^), see main text). In dashed line, we plot the expected 
numbers of modules in the natural partitions. In the lower figure, we 
plot the measures of robustness (V)ai go (t), (V)net(t) and (V)qf(£). 
Natural partitions are robust and associated to low values of (V), 
while in-between values of t are characterized by peaks of (V) 
(except for (V%i go (t) at me jump from 4 modules to 1 module). 
Vertical lines indicate peaks in (V) and the behavior of the system 
at t — 1 (modularity). 



B. College football 

This real- world network is made of 115 football teams 
that are connected if they have played a regular-season 
game [4]. Because games are more frequent between 
members of the same conference than between members 
of different conferences, one expects a natural partition 
into 12 communities, corresponding to the 12 confer- 
ences of the championship. For this network, we use 
the parameters K = 100, T = 100 and A = 5. One 
observes (see Fig. 3) a clear plateau where partitions are 




time 



Fig. 3. Analysis of the college football network. In the upper figure, 
we plot the average number Nt of modules and its standard deviation 
a t when performing 100 optimizations with random orderings. In 
dashed line, we plot the number of modules in the expected partition 
(12 communities). In the lower figure, we plot the measures of ro- 
bustness (V) a \g (t), (V)net(t) and (V)qf(£)- The only time window 
where (V) vanishes corresponds to a partition into 12 communities. 
Interestingly, t = 1 (indicated by a vertical line) is not particularly 
robust, i.e., modularity does not uncover modules at the appropriate 
scale. As shown in the upper figure, for each value of t, the numbers 
of modules in the 100 optimal partitions are close to their average 
Nt (small values of a t ) even when the partitions are not robust. 



made of 12 communities and where the three versions of 
robustness are vanishingly small. It is interesting to note 
that t = 1 does not belong to this plateau. Modularity 
optimization thus fails to uncover a robust partition and 
provides an inappropriate representation of the system. 
This intrinsic problem of modularity has already been 
observed in benchmarks |34l and empirical networks 

ma. 

VI. Discussions 

In this article, we have focused on the detection of 
non-overlapping modules in multi-scale networks. These 



networks are made of different levels of organization 
and are typically (but not necessarily) hierarchical, in 
the sense that the system is made of modules, which 
themselves are made of sub-modules, etc. We have 
shown that modularity optimization is not a satisfac- 
tory method to uncover modules in general, because 
modularity optimization reveals communities at scales 
that are not automatically compatible with the system 
organization. It is therefore necessary to incorporate a 
resolution parameter to modularity in order to adjust the 
characteristic size of the modules and to uncover the true 
modular organization of a network. Three different multi- 
resolution quality functions Q 7 , Q r and Q t have been 
presented. They are all equivalent up to a linear transfor- 
mation and include modularity as a particular case when 
the resolution parameters are 7 = 1, r = and t = 1. 
It is important to keep in mind that multi-resolution 
quality functions have the same limitations as modularity 
when the resolution parameter is fixed [ 35 1 and that the 
possibility to tune this parameter is essential. No value 
of the resolution parameter is a priori better than another 
one and additional tests are therefore needed to uncover 
significant scales of description. Contrary to what is 
usually believed, modularity is thus an ordinary instance 
in the set of multi-resolution quality functions and there 
is no deductive reason to prefer it. Important values 
of the resolution parameter have instead to be selected 
by considering the robustness of the detected partitions. 
Our analysis suggests to develop proper statistical tests 
and to combine the information obtained from different 
measures of robustness in order to better comprehend 
the modular organization of complex networks. 
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